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Abstract: We calculate the hadronic matrix elements to 0(p 4 ) in the chiral expansion 
for the (AS = 1) K° -* 2vr decays and the (AS = 2) K°-K° oscillation. This is done 
within the framework of the chiral quark model. The chiral coefficients thus determined 
depend on the values of the quark and gluon condensates and the constituent quark 
mass. We show that it is possible to fit the AI = 1/2 rule of kaon decays with values 
of the condensates close to those obtained by QCD sum rules. The renormalization 
invariant amplitudes are obtained by matching the hadronic matrix elements and their 
chiral corrections to the short- distance NLO Wilson coefficients. For the same input 
values, we study the parameter B^ of kaon oscillation and find = 1.1 ± 0.2. As an 
independent check, we determine Bx from the experimental value of the K^-Ks mass 
difference by using our calculation of the long-distance contributions. The destructive 
interplay between the short- and long-distance amplitudes yields Bk = 1.2 ± 0.1, in 
agreement with the direct calculation. 
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1 Introduction 

The physics of kaons is an important testing ground for our understanding of low-energy 
QCD. The hadronic matrix elements for the nonleptonic decays must satisfy the non- 
trivial constraint of the AI = 1/2 rule that shows the striking enhancement of the decay 
width in which the two final pions combine in the isospin 1 = state with respect 
to the 1 = 2 state. The interplay between the renormalizat ion-group evolution of the 
Wilson coefficients and the matrix elements of the relevant quark operators in both K°-K° 
mixing and K° — > 2n decays determines the theoretical predictions of the CP- violating 
parameters e and e'/e, the knowledge of which is essential to the standard model. 
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The recent progress in the next-to-leading-order (NLO) computation of the Wilson 
coefficients of the AS = 1 [1] and AS = 2 [2] effective lagrangian at the quark level 
makes more urgent the need to bring under better control the non-perturbative estimate 
of the corresponding hadronic matrix elements. Different approaches have been pursued, 
ranging from lattice simulations [3], 1/A^-expansion [4], phenomenological approaches [5] 
dispersion-integral techniques [6] and low-energy QCD modeling [7, 8] . A particular case 
of the latter, the chiral quark model (xQM) [9, 10] has been analyzed in detail in a 
series of papers [11, 12, 13] with the encouraging result of fitting the AI = 1/2 rule by 
a consistent and reasonable choice of the input parameters. These are essentially three: 
the quark condensate, the gluon condensate and M — a parameter of the model which 
corresponds to a typical constituent quark mass. 

In a previous paper [11], we have calculated the coefficients of all the 0(p 2 ) terms of 
the AS = 1 chiral lagrangian and included the corresponding chiral loop renormalizations 
(which are of 0(p 4 )). The present paper completes our previous analyses 

• by including in the hadronic matrix elements the complete NLO 0(p 4 ) corrections; 

• by updating the short-distance calculation of the Wilson coefficients according to 
the most recent determinations of a s and m t . 

The determination of the K — > 2tc matrix elements at 0(p 4 ) is a non-trivial compu- 
tation. To our knowledge it is the first time that the whole of these hadronic matrix 
elements are estimated to this order by any technique at all (for previous discussion, see 
refs. [14, 15]). 

The introduction of these corrections allows us to determine a consistent range of 
values of the input parameters for which the AI = 1/2 rule is reproduced with an accuracy 
at the 20% level. We find the best fit to this order in the chiral expansion for values of 
the condensates close to those derived by QCD sum rules and a value of M in agreement 
with independent estimates based on radiative kaon decays [16]. A coherent picture is 
thus provided: all hadronic matrix elements are computable for a common set of input 
parameters and no ad-hoc assumption is necessary to fit the AI = 1/2 rule. 

The hadronic matrix elements thus found are eventually matched together with one- 
loop chiral corrections to the corresponding NLO Wilson coefficients. The matching scale 
is a delicate choice in so far as it can neither be too high, in order for chiral perturbation 
theory to remain valid, nor too low, for the Wilson coefficients to be reliable. We choose 
a scale of 0.8 GeV that we identify with the chiral symmetry breaking scale. Even though 
at this scale a s is already rather large, we have verified that the computed observables 
change by no more than 30% when moving from LO to NLO order Wilson coefficients. 
By considering the shift in the value of a s at /i = 0.8 GeV when going from the NLO 
to the next order, we can estimate that the neglect of next-to-NLO corrections in the 
running of the Wilson coefficients may affect our predictions at the 10% level, which is 
within the "systematic" error we assign to our analysis. 

Once the input parameters of the model have been fixed, other quantities of kaon 
physics can also be estimated. In particular we discuss the deviation from the vacuum 
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insertion approximation (VSA) in K°-K° oscillations that is parameterized by the scale 
independent parameter Bk- The inclusion of the 0(p 4 ) corrections to the AS = 2 chiral 
lagrangian provides us — for the input parameters identified by the AI = 1/2 selection 
rule — with a central value of Bk — 1.1. 

Contrarily to what one might think, such a large value is not in conflict with the 
experimental determination of the K^-Kg mass difference AM^g. To prove this point, 
we have computed the long-distance mesonic contributions that arise from the double 
insertion of the AS = 1 chiral lagrangian, and find that they are about 20-30% of 
the whole mass difference and of the opposite sign with respect to the short-distance 
part [17]. By requiring that the AM LS thus calculated fits the experimental value, we 
obtain Bk ~ 1.2 in striking agreement with the direct calculation. 

We thus present two independent estimates of B K within the xQM approach: 

• Bk = 1.1 ± 0.2 from direct calculation, 

• B K = 1-2 ± 0.1 from AM LS) including long-distance contributions, 

where the errors include a flat variation of all relevant input parameters. According to 
these results, values of B K smaller than one are disfavored in the xQM- 

These results can then be applied to the prediction of the direct CP violating param- 
eter e'/e . Because of the importance and the complexity of such a calculation, we have 
decided to present it in an independent paper [18]. 

In order for the present paper to be as self-contained as possible, we have included 
in the following two subsections the relevant lagrangians and a brief introduction to the 
xQM. Such introductions summarize those of our previous papers, to which we refer the 
reader for further details and references. 



1.1 The Quark Effective Lagrangian 

Let us introduce our notation by recalling that the AS = 1 quark effective lagrangian at 
a scale jj, < m c can be written as [19] 



The Qi are effective four-quark operators obtained by integrating out in the standard 
model the vector bosons and the heavy quarks t, b and c. A convenient and by now 
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standard basis includes the following twelve quark operators: 



Qi = (s a u ( 3) y _ A (upd a ) Y _ A , 

Q2 = {su)v- A (ud) y _ A , 

Q3,5 = (sd) y _ A Eg (^)v T A > 

Qifi = (~Sadf3) Y _ A Y,q(q/3 ( la)vTA, 

Q 7 ,9 = § (sd) v _ A Eg e q (qq) Y±A , 

<58,10 = ^(s a d^) Y _ A Y^q^q(Q^Qa)v±A, 

Q11 = 1^2 s [m d (l + 75) +m s (1 - 75)] a ■ Gd, 

Q12 = 1^2 s [m d (1 + 75) + m s (1 - 75)] a-Fd, 



(1.2) 



where a, /3 denote color indices (a, (3 — 1, . . . , iV c ) and e g are quark charges. Color indices 
for the color singlet operators are omitted. The subscripts (V ±A) refer to 7 /i (l±7s). We 
recall that Qi,2 stand for the ^-induced current-current operators, Q3-6 for the QCD 
penguin operators and Q7-10 f° r the electroweak penguin (and box) ones. The quark 
operators Qn,n, involving the gluon and photon fields, are the dipole penguin operators 
whose K° — > 7T7T matrix elements arise at 0(p A ) [14]. 

Even though not all the operators in eq. (1.2) are independent, this basis is of particular 
interest for the present numerical analysis because it is that employed for the calculation 
of the Wilson coefficients to the NLO in a s [1]. 

In the present paper we will mainly discuss the features related to the first six operators 
in (1.2) because the electroweak penguins Q7_io have their contributions suppressed by the 
smallness of their CP conserving Wilson coefficients, while Q11-12 have very small matrix 
elements [14]. The electroweak operators play however a crucial role in the discussion of 



The functions Zi(ji) and yi(fi) are the Wilson coefficients and the Kobayashi-Mas- 
kawa (KM) matrix elements; r = —V t dV t * s /V u dV* s . The numerical values of the Wilson 
coefficients at a given scale depend on a s . We take the most recent world average [20] 



We match the Wilson coefficients at the m w scale with the full electroweak theory by 
using the LO MS running top mass m t (mw) = 177±7 GeV which corresponds to the pole 
mass m? t ole = 175 ±6 GeV [21]. For the remaining quark thresholds we take mb(mb) = 4.4 
GeV and m c (m c ) = 1.4 GeV. 

Similarly, the effective AS = 2 quark lagrangian at scales ji < m c can be written as 

£as=2 = -C S 2(tj) Qs2{jj) 



e'/e [18]. 



a s (m z ) = 0.1189 ±0.0020, 
which at the NLO corresponds approximately to 

AL 4 L = 340 ± 40 MeV . 



(1.3) 



(1.4) 



4tt 2 



\ 2 c r)iS{x c ) + \ 2 t rj 2 S{x t ) + 2X c X t rj 3 S(x c ,x t ) b(fi)Q S 2(n) (1.5) 
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where Xj = Vj d V* s and Xi = mf/m^. We denote by Q S 2 the AS = 2 local four quark 
operator 

Qs2 = {s L Yd L ){s Ll ,d L ) ) (1.6) 

which is the only numerically relevant operator in this case. 

The integration of the electroweak loops leads to the functions S(x) and S(x c , xt), the 
exact expressions of which can be found in ref. [22]. They describe the dependence of the 
AS = 2 transition amplitude on the masses of the charm and top quarks in the absence 
of strong interactions. 

The short-distance QCD corrections are encoded in the coefficients rji, r] 2 and i] 3 and 
b(fx). The r]i coefficients represent the renormalization effects down to the scale m c . 
They are functions of the heavy quarks masses and of Aq CD . These QCD corrections are 
available to NLO in the strong coupling [2]. The scale-dependent function describes 
the overall running below the charm threshold and it is given by 

%) = [a. (^)]- 2/9 (l - J 3 ^) , (1.7) 

where J3 depends on the 75-scheme used in the regularization. The naive dimensional 
regularization (NDR) and the 't Hooft-Veltman (HV) scheme yield, respectively: 

■h ~ 162 and J 3 - ig2 . (1.8) 
1.2 The Chiral Quark Model 

In order to evaluate the bosonization of the quark operators in eqs. (1.2) and (1.6) 
we exploit the xQM approach which provides an effective link between QCD and chiral 
perturbation theory. 

The xQM can be thought of as the mean field approximation of the extended Nambu- 
Jona-Lasinio (ENJL) model for low-energy QCD. A detailed discussion of the ENJL model 
and its relationship with QCD — as well as with the xQM — can be found, for instance, 
in ref. [8]. 

In the xQM, the light (constituent) quarks are coupled to the Goldstone mesons by 
the term 

4^ M = —M (q R Sg L + q L tfq R ) , (1.9) 
where q T = (u, d, s) is the quark flavor triplet, and the 3x3 matrix 

£ = exp^yn(:r)^ (1.10) 

contains the pseudoscalar octet Tl(x) = J2 a X a ir a (x)/2, (a = 1,...,8). The scale / is 
identified at the tree level with the pion decay constant f w (and equal to fx before chiral 
loops and higher order corrections are introduced). In the xQM, the pion decay constant 
is given by a logarithmic divergent quark loop integral , which is numerically identified 
with /, such that /+ EE /f// = 1. 
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The xQM has been discussed in several works over the years [9, 10]. We opted for the 
somewhat more restrictive definition suggested in ref. [10] (and there referred to as the 
QCD effective action model) in which the meson degrees of freedom do not propagate in 
the original lagrangian. 

The QCD gluonic fields are considered as integrated out down to the chiral breaking 
scale A x , here acting as an infrared cut-off. The effect of the remaining low- frequency 
modes are assumed to be well-represented by gluonic vacuum condensates, the leading 
contribution coming from 

(-GG). (1.11) 

7T 

The constituent quarks are taken to be propagating in the fixed background of the soft 
gluons. This defines an effective QCD lagrangian £^ D (A X ), whose propagating fields are 
the u, d, s quarks. The strong xQM lagrangian is therefore given by 

£ X QM = £qcd( A x) + ^xQM • I 1 - 12 ) 

The £ x qm interpolates between the chiral breaking scale A x and M (the constituent 
quark mass). The three light quarks (u, d, s) are the only dynamical degrees of freedom 
present within this range. The total quark masses are given by M + m q where m q is the 
current quark mass in the QCD lagrangian. The Goldstone bosons and the soft QCD 
gluons are taken in our approach as external fields. A kinetic term for the mesons, as 
well as the complete chiral lagrangian, is generated and determined by integrating out 
the constituent quark degrees of freedom of the model. By combining eq. (1.12) with 
eqs. (1.1) and (1.6) one may obtain the AS = 1 and AS* = 2 weak chiral lagrangians 
as effective theories of the xQM. In the matching process, the many coefficients of the 
chiral lagrangian are determined — to the order 0(a s N c ) in our computation — in terms of 
M, the quark and gluon condensates. We neglect heavier scalar, vector and axial meson 
multiplets. 

In conventional chiral perturbation theory the scale dependence of meson loops renor- 
malization is canceled by construction by the 0(p 4 ) counterterms in the chiral lagrangian. 
While in our approach this is mantained for the strong sector of the chiral lagrangian, 
the tree-level counterterms of the weak sector are taken to be \x independent and a scale 
dependence is introduced in the hadronic matrix elements via the meson loops, evaluated 
in dimensional regularization with minimal subtraction. This scale dependence is eventu- 
ally matched with that of the Wilson coefficients. Tree level counterterms acquire a scale 
dependence at the next order in the chiral expansion, via meson loop renormalization. 

A more detailed discussion about our approach including the complete determination 
of the AS = 1 chiral lagrangian at 0(p 2 ) can be found in refs. [11, 14]. 

2 Building Blocks 

Within the xQM, matrix elements of the quark operators Qi-io can be calculated in the 
factorizable approximation as products of two currents or two densities. Non-factorizable 
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matrix-elements are proportional to the gluon condensate. In this way the hadronic 
matrix elements of the relevant operators are constructed in terms of products of two 
building blocks. These elementary blocks must be computed to the appropriated order 
in p 2 and m q to yield the 0(p A ) matrix elements. 

By integrating over quark loops we obtain the following matrix elements of quark 
densities, based on the lagrangian in eq. (1.12): 



s^u\K + (k)) = iy/2 



(n+( P+ )\sd\K + (k)) = 



m -k^(f ++ k2 



f 



2M 



A 2 



+ (m s + m u ) f /+ + 3 



+ 



f_ 

M 

(qq) 



M 2 



P 



(m 2 s +m 2 u ) y+- Ql j^rJ - A 



2M' 



2A 2 



16MA 2 



1 (P< + 2py + 5q <) 



(m 8 + m d ) 



"2 Vtly 
1 

~M 



U + 



3P 2 + 5q 2 



(m s - m d ) 



q-P 

A 2 

x 



, 2 2N fr M ^ 

[m s + m s m d + m d ) I /+ - 6— 

M 2 

- 6 m u (m s + m d + m u ) — 



(2.1) 



(2.2) 



where (qq) is the quark condensate for zero current quark mass. In the calculation of 
the four-fermion matrix elements at 0(p 4 ) also the density amplitude (7T7t| s'-f^d \K°) con- 
tributes. A discussion on its evaluation is given in appendix A. 
The corresponding quark current matrix elements are given by: 



^f(l- l5 )u\K + (k)) 



(n+(p + )\s^(l- l5 )d\K + (k)) 



+ 



m s + m u 



2M 



U ~ 12 



M 2 

A 2 

x 



u + 



P 2 + 3g 2 
2A 2 



q ■ P 

"K 2 

X 



-q 



3M(m s + m d + 2m u ) 



m s - m d 
2M 



M 2 
x J 



(2.3) 



(2.4) 
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(0\s^T a (l- l5 )u\K + (k)) 
(n + ( P+ )\sYT a (l- l5 )d\K + (k)) 



:V2 



16tt 2 / 



16tt 2 / 2 



Gl T A^{k) 



(2.5) 
(2.6) 



where q = k — p + and P = k + p + , k being the incoming momentum of the kaon field 
and p + the outgoing momentum of the pion field. The chiral symmetry breaking scale is 
taken as A x = 2n^j6/N c f ~ 0.8 GeV, while /+ = ffl/f = 1 is the vector form factor at 
zero momentum transfer. Finally, G a VT with a = 1,...,8 is the usual SU(3) C gluon field 
tensor and T a are the SU(3) C generators normalized as Tr T a T b = 5 ab /2. 

The eqs. (2.5)-(2.6) represent the gluonic corrections, which are computed by using 
one-gluon dressed quark propagators and color Fierz transformations on the four-quark 
operators (the SU(3) C index a is to be summed over in the full matrix element). They 
contribute to the non-factorizable part of the hadronic matrix elements via the relation 



which allows us to express our results in terms of the gluonic vacuum condensate. 
The two gluonic form factors are given by: 

k 2 



(2.7) 



1 + 



TO,, 



and 



e(P, n, v, r) 



+ i- 



1 + 



.TOo 



TO 



4M 



6M 2 2M 



(2.8) 



P 2 



+ 



12M 2 AM 2 



-e{q, //, u, t)—- - e(P, q, fi, v)-r—^ + e{P, q, fi, r 



6M 2 



12M 2 



12M 2 



2 P 



120M 2 

TO 



29P 2 llo 2 
+ 



240M 2 80M 2 



9 HT Qr9piV 



12M 

m d 

12M 

m u 

12M 



e [(5P + q), /x, u, r] + - \{P V - 19g„) - (v -> r)] 
{e [(5P - g), i/, r] - | [{P v + 19g v ) <fo r - (u - r)] 
|l2e (P, /x, !/, r) - ^ [14^^ T - (i/ -> r)] } 



(2.9) 



where e(k,n,is,r) = k a e aflUT . 

The procedure by which it is possible to determine the hadronic matrix elements and 
the corresponding chiral coefficient to the desired order for all the relevant operators by 
means of the given building blocks is discussed in refs. [11, 14], to which we refer the 
interested reader. 
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Before proceeding it is important to identify what parts of the building blocks are 
absorbed by the renormalization of the parameters and fields of the lagrangian. To this 
we now turn. 



3 Wave-Function and Coupling- Constant Renormalizations 

To the order 0(p 4 ) at which we are working, both the wave-function and the meson 
decay constant renormalizations must be included. They give sizeable contributions that 
cannot be neglected. In our computation, the renormalization comes in two parts which 
are conceptually distinct. On the one hand we have the chiral renormalization generated 
by the meson loops, to which we have to add the renormalization that is specific of the 
xQM and originates from the expansion of the building blocks beyond the leading order. 
In refs. [11, 12, 13] only the chiral loop renormalizations were included. 



3.1 Wave-Function Renormalizations 

The wave-function renormalizations which arise in the xQM from direct calculation of 
the K — > K and n — > n propagators are given at 0{p 2 ) by: 

Z " = 1 - 2 AJ +6 AJ 

Z. = l-^ + 6 M(m ;+ m "» (3.1) 

The complete 0(p A ) expressions are given in appendix B. The renormalizations above are 
added to those induced by the one-loop chiral corrections. 



3.2 Renormalization of f K and f n in the xQM 



The building block in eq. (2.3) gives directly the corrected (unrenormalized) version of 
the coupling constant / yelding (for instance for f K ) 



fl = f 



i _|_ H^k + m s + m t 



A 2 

x 



2M 



1 - 12 



M 2 

A 2 " 

x 



(3.2) 



However, it is only after inclusion of the wave-function renormalization in eq. (3.1) that 
we find the complete xQM expression: 



fK = \/Z K f K = f 



u 



1 + 



2M 



1-6 



(3.3) 



and similarly for f n . Notice that f^ K coincide at 0(p 4 ) with the physical decay constants 
f n> K only after inclusion of the proper chiral loop renormalizations (see next subsection). 

The consistency of the entire procedure can be verified by considering the renormal- 
ization of the charged form factor f+{q 2 ) in the coupling of mesons to an external vector 
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field (photon or vector boson). This form factor can be read off eq. (2.4). After renor- 
malization we find that the on-shell form factor is 



X 

which correctly preserves the current-conservation condition /+(0) = 1. 

By matching the expressions above for the renormalized fx and f n with the chiral 
lagrangian results we obtain the xQM determinations (see also appendix B) 

L 4 = , (3.5) 
f 4 / M 2 \ 

L - = -mm I 1 - 6 a*)- (3 - 6) 

where by PCAC 

fe)= _^L = _^L. (3.7) 
m s + m u m d + m u 

The above equations are obtained at the constituent quark mass scale M, where the 
quark fields are integrated out and the standard chiral lagrangian arise as the effective 
theory of the xQM. The parameter L 4 is vanishing in the xQM up to higher order gluon 
condensate contributions. The usual numerical determinations of the renormalized L 4 
and L 5 in chiral perturbation theory must then be run down to this scale in order to be 
compared to the result of the xQM, after taking into account the different subtraction used 
in the renormalization of the chiral loops. We address this issue in the next subsection. 

3.3 The Chiral Coupling / at the One-Loop Order 

The one-loop expressions for f n and f K allow us to determine the value of the chiral 
coupling / to be used in the NLO calculation. In order to compare eq. (3.3), and the 
analogous one for f n , with the physical values of the decay constants, the 0(p 4 ) renormal- 
ization induced by chiral loops has to be included. In terms of the tree-level counterterms 
of the 0{p A ) strong chiral lagrangian, one obtains the equations [23] 

f - f w f - 2h(m n , fj) - h(m K , fi) + AmlLM + (4ml + 8m 2 K )L 4 (fi) = , (3.8) 

for the pion decay constant f n , and 

3 3 3 

f ~ fxf - -^h(m w , fi) - -h(m K , /i) - -h(m v , fi) + Am 2 K L 5 (/j,) + (Ami + 8m 2 K )L 4 (/j,) = , 

(3.9) 

for fx, where, by keeping only the logarithmic part from the loop integrals, 

2 2 
TTl 771 

H m > A*) = 32^5 ln ■ ( 3 - 10 ) 

This prescription defines the renormalized couplings L 5 and L 4 . Once L 4 is given, these 
equations may be solved for L 5 and /. 
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The entire procedure contains many uncertainties, the input value of L 4 being one of 
them. According to the original derivation of ref. [23] we take L^(m v ) = 0.0 ± 0.5. Using 
the known anomalous dimensions of the couplings L; L [23], we run L 4 up to the matching 
scale \x = A x ~ 0.8 GeV where we find 

L 4 (A X ) = -(0.3 ±0.5) x 1(T 3 . (3.11) 

For this value, the solution of the system of eq. (3.8) and eq. (3.9) at the scale /i — A x 
gives 

L 5 (A X ) = (1.2 ± 0.4) x 10" 3 and / = 0.086 ± 0.013 GeV , (3.12) 

where the errors come from uncertainties in the input parameters, mainly those for L 4 . 
By solving the same equations at the scale M ~ 0.2 GeV, where 

L 4 (M) = (0.8 ± 0.5) x 10~ 3 , (3.13) 

we obtain 

L 5 (M) = (4.3 ± 0.4) x 10~ 3 (3.14) 

and / = 0.086 ± 0.013 GeV, as expected from the scale-independence of /. 

A consistent solution of eqs. (3.8)-(3.9) is also obtained by using the MS result for 
h(m, /i), that is 

=^( b, S (3 - i5) 

In this respect, it is important to remark that the determination of L5 and L 4 depend 
on the subtraction scheme employed in the renormalization of the chiral loops (while 
the value of / is by construction independent on it). Only by consistently adopting the 
same renormalization prescription the physical quantities remain invariant. By taking 
/ = 0.086 ± 0.013 GeV we obtain the following ITS results: 

Z 4 (M) = (0.4 ^) x KT 3 , (3.16) 

L 5 (M) = (3.2 ±0.3) x 10~ 3 , (3.17) 

and 

Z 4 (A X ) = (-0.71™) xlO- 3 , (3.18) 

L 5 (A X ) = (0.1 ±0.3) x 10^ 3 . (3.19) 

The value of L A in eq. (3.16) is consistent with the xQM determination in eq. (3.5). The 
matching of the expression in eq. (3.6) with the result of eq. (3.17) gives a value of (qq) 
at pi — M of (—185 ± 10 MeV) 3 . The larger values of the quark condensate that can be 
extracted at the scale A x from eq. (3.19) are consistent with the range of values of (qq) 
that we obtain at /i = A x from the fit of the AI = 1/2 rule, computed in the MS scheme. 

Even though the actual values for the parameters (qq) and (a s GG/n) that fit the 
selection rule will vary depending on the choice of / in the range of eq. (3.12), we have 
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verified that the predictions of the other observables remain rather stable when varying 
the input parameters accordingly. In particular, the value of Bx varies by a few percents 
by changing / within the range of eq. (3.12), provided we use the corresponding values of 
the input parameters that fit the AI = 1/2 rule. We will therefore take the central value 
of / as a fixed input in all numerical estimates. 



4 Hadronic Matrix Elements 

The chiral lagrangian coefficients to order 0{p 2 ) were computed in previous papers: 
refs. [11] and [17] for AS = 1 and AS = 2 transitions respectively. Here we report 
the NLO corrections to the matrix elements (Qi)f LO which are obtained by means of the 
building blocks in section 2. We define m u = rrid = m, the average of the u and d current 
quark masses, consistently to m^i = m^o = m n . 

There is no mass renormalization to this order because all masses in the matrix ele- 
ments of the operators Qi-e originate from the external momenta which are defined to 
be on the physical mass-shell. As expected, constant terms in the densities, which enter 
in the determination of the operators (5s,6> cancel in the complete matrix elements. 

The total matrix elements have the form 



where Qi are the operators in eq. (1.2), 

(Qi)! = ((n^^QilK 



NLO 



(4.1) 



(4.2) 



After including the wavefuntion renormalization, the matrix elements are expanded to 
0(p 4 ), discarding higher order terms. The functions af(/z) represent the scale depen- 
dent meson-loop corrections, including the mesonic wavefunction renormalization. They 
are defined as the isospin projections of the af~(/i) and a°°(/i) corrections computed in 
ref. [11], properly rescaled by factors of f w / f in order to replace /„.—>/ in the NLO eval- 
uation. The scale dependence of the NLO part of the matrix elements is a consequence of 
the perturbative scale dependence of the current quark masses which enter at this order. 

The LO matrix elements (Qi)j , ° can be found in ref. [11], all occurences of f n beeing 
replaced by /. The NLO 1 = and 2 contributions to the K° — > tttt matrix elements are 
given by: 



NLO 



(Qi)? LO 



(Q^ LO 
(Q^ LO 



3 

V2 



X 



- l -x 

3 



(Qi)? LO , 



- 1 + !>-atW g 



N c 



-2 + ^-^ G + 3 7G ) 



(4.3) 
(4.4) 

(4.5) 
(4.6) 
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(Qs)Z LO 
(QX LO 

m N o LO 

(QX LO 



■^-X [(3 - 5 {GG) {(3 G - 7 G )] , 

(QX LO - (QX LO + (Qs)Z LO , 



where 



X = V3f (m\ - ml) 
and the gluon-condensate correction 5^cg) is given by 

N c (a s GG/n) 



>(GG) 



2 16tt 2 / 4 

The quantities (3,0g,1 an d 7g are dimensionless functions of the mass parameters: 



(4.7) 
(4.8) 
(4.9) 

(4.10) 
(4.11) 

(4.12) 
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7 
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+ Sfl-12^ 
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6M 2 
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12M [m\ - ml) M 
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m\ m 2 K {2m s + 6m) — m\{m s + 7m) 



2MA 2 X 



(m 2 K - ml) A 2 



+ 



2m(m s — m) 
M(m 2 K - ml) 

2 



2f + -6 



M 2 



+/ 



m^ + m 2 ^ (m s + 5m) ^ 2m(m s — m) 



AM 2 



2M 



rn 



K 



mi 



A 2 

x 



m\ m 2 K {m s + 3m) — 2m\{m s + m) 



m 



7g 



A 



2M 

m, — m 



(m|- - m 2 ) 



A' 
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Similarly, for the AS = 2 matrix element {K°\Q 2S \K°) we write 



{QM) = Z K [(Q 2S ) LO + (Q 2S (ri) NLO }+a 2S (fi) 



with 



and 



(Q2s) L ° = m\f 



2 e2 



NLO 



m 2 K f { (l + 



1 + Tr(l- <W 



N, 



m 



A 



+ 



m s +m 



M 



1 - 12 



A 2 " 

X 



${gg) ( m\ 
N r \3M 2 



m s + m 
M 



(4.13) 
(4.14) 



(4.15) 
(4.16) 

(4.17) 
(4.18) 



(4.19) 
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where the scale dependence enters through the perturbative running of the current quark 
masses. 

The chiral-loop correction to (K°\Q 2S \K°) , including the meson wave-function renor- 
malization, is given by 



1 + ^ 



Hgg), 



a 2S {n) = m 2 K f 2 K 

where, having normalized the 0{p A ) amplitude to iri^f^, 

r 2S (fj,) = 1.45 + 0.676 ln/i 2 



(4.20) 



(4.21) 



and \i is in units of GeV. The given numbers depend only on the octet meson masses and 
are obtained for the values given in appendix C. 

By subtracting the MS chiral loop renormalization of f 2 — > /J- at the scale fj, one 
obtains [17] 

r 2S {n) = 0.728 + 0.373 ln/i 2 , (4.22) 

which shows that a large part of the meson-loop correction goes into the renormalization 
of/*. 

Everywhere, in the NLO as well as in the leading order matrix elements, the coupling 
constant / must be understood as given by the values discussed in section 3.3. This 
replacement in the NLO terms is optional, given that the difference is of 0(p 6 ). For 
convenience, in our numerical analysis we replace all occurrences of / with its one-loop 
value. 

Another remark concerns the quark condensate which appears in the building blocks of 
the quark densities. While all factorizable gluonic corrections are absorbed in the physical 
definition of the chiral coupling and the quark condensate, there are mass dependent 
contributions to the latter that we have not included in eq. (4.15). At the nedeed 0(m q ), 
the relation between the condensate (qq) that we use as a parameter in our analysis, and 
the condensate (qq) mq calculated at non-zero current quark masses is given by 



m 



(qq) ma = (qq) + -jf (qq) + 2Mf 2 f. 



(4.23) 



The correction is negligible for q = u,d and remains small also for (ss) due to a partial 
cancellation between the two terms in square brackets. In what follows we will always 
refer to the flavor independent parameter (qq). 



4.1 Chiral Loop Corrections 

The one- loop chiral corrections ae(fi) are included at this stage operator by operator on 
top of the xQM estimate of the matrix elements, as discussed in refs. [11, 17]. Let us 
recall that throughout our analysis we use dimensional regularization and the modified 
minimal subtraction scheme as it is done for the Wilson coefficients. For this reason the 
numerical values for the AS = 1 chiral loop corrections quoted in ref. [11] differ from those 
of ref. [15] where a non-minimal subtraction scheme is employed. The two computations 
are however in agreement once the scheme dependence is taken into account. 
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5 The AI = 1/2 Selection Rule 



As discussed in the introduction, in order to restrict the possible values of the input 
parmeters M, (qq) and {a s GG/n) we study the AI = 1/2 selection rule which characta- 
rizes the 1 = and 1 = 2 amplitudes of the non-leptonic kaon decays. 

Is is convenient to write the generical amplitude for a kaon to decay into a two pion 
final state of isospin I = 1,2 as 

Ai(K -> tttt) = Ai exp i (5/) (5.1) 

where the phase 5i comes from the final state interactions in the channel I. From eq. (1.1) 
and that fact that Imr < 1 we can write the amplitude Ai as 

cos 5/ 

Experimentally the phases 5i are obtained in terms of the ix-it S-wave scattering 
lenght [24] at the m K scale. The values so derived give to a few degrees uncertainty 
5q ~ 37° and 5 2 ^ —7°, thus obtaining with good accuracy 

cos5 — 0.8 

cos<5 2 ~ 1.0. (5.3) 

As a consequence the 1 = amplitude includes a 20% enhancement from the rescattering 
phase. 

The existence of 1m Aj ^ signal the possible presence of direct CP violation in 
the K — > 7T7T decays, while uj = |^4 2 |/|-^o| — ReA 2 /ReA = 1/22.2 represents the 
AI =1/2 selection rule known experimentally for more than forty years [25]. According 
to our conventions for the isospin amplitudes we have ReA = 3.3 x 10~ 7 GeV and 
ReA 2 = 1.5 x 10~ 8 GeV. Using the notation of eq. (1.1) we can write the CP conserving 
amplitudes as 



ReA2 = ^zY^l^Zi Re(Q l ) 2 (1 + *W) » ( 5 - 4 ) 
V2 coso 2 i 

where Q v +r/ = 0.25 ±0.05 represents the isospin breaking effect due to the 7i -r]-r]' mixing. 
For notational convenience henceforth we identify Re Aj with Aj. 

For a long time the explanation of the AI = 1/2 rule has been a mystery. It was 
found some twenty years ago that perturbative QCD qualitatively worked in the right 
direction [26, 27] — even though the effect was too small. More recently, it has been 
shown shown that non-factorizable contributions and chiral loops also work in the right 
direction [4, 7, 15]. In a previous paper [12], we argued that the selection rule could be 
accomodated within the context of the xQM to 0(p 2 ) with the inclusion of the meson 
loops. Here we present the complete 0{p A ) calculation. 
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To this order in the chiral expansion there are several constraints to be satisfied. In 
particular, the stability of the solution as well as its 75-independence restrict the possible 
values of M. From Figs. 1 and 2, obtained for the central values of the parameters given 
in appendix C, we see that in order to keep the numerical results for Aq and A2 within 
20% from their experimental values we must constrain M in the range 

M = 180^220MeV (5.5) 

The range thus identified is in good agreement with that found by fitting radiative kaon 
decays [16]. We notice that the NLO corrections improve the 75-stability, especially for 
A 2 . Henceforth we will use for the numerical discussion the HV results. 




0.16 0.18 0.2 0.22 0.24 

M (GeV) 



Figure 1: Aq(K° — > -kit) as a function of M in the HV (gray) and NDR (black) schemes 
respectively. The horizontal line indicates the experimental value. 



2.5 




0.16 0.18 0.2 0.22 0.24 

M (GeV) 



Figure 2: A2(K° — ► tttt) as a function of M in the HV (gray) and NDR (black) schemes 
respectively. The horizontal line indicates the experimental value. 

Given the above range for M, we can study the dependence of the selection rule on 
the other xQM parameters. We find that for central values of M and Aqq D , namely 
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Figure 3: Aq and A2 in units of GeV as functions of {a s GG/ir) (vertical spread) and {qq) 
(horizontal spread) for fixed M = 200 MeV and central value of Agjj D . The cross-hairs indicate 
the experimental values. Black and gray lines correspond to m s (0.8 GeV) = 240 MeV and 200 
MeV respectively. The ranges of the gluon and quark condensates are those discussed in the 
text. The slight slope in the curves is due to the {qq) dependence in the electroweak operators. 



M = 200 MeV and A { q ] cd = 340 MeV, the AI = 1/2 rule is reproduced with a ±10% 
accuracy provided 

{a s GG/n) = (334 ± 4 MeV) 4 , 

and 



(qq) 



240 ±w M eV 



(5.6) 
(5.7) 



as it is shown in Fig. 3. 
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Figure 4: Aq and A2 in units of GeV as functions of {a s GG/ir) (vertical spread) and {qq) 

220 MeV and the extreme 
300 MeV) 



(horizontal spread) in the ranges of eqs. (5.6)-(5.7) for m s (A x ) 
values of M and pSq CD . Black and gray lines correspond to M = 205 MeV {Aq CD 



and 197 MeV (A^ D = 380 MeV) respectively. 

Given the results above, we include the dependence on Aq C d- Requiring that the rule 
is finally reproduced with a ±20% approximation, for A^q CD in the range of eq. (1.4) we 
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obtain a sharp constraint on M (Fig. 4): 

M = 200 ±1 MeV . (5.8) 

In order to exhibit the impact of the NLO corrections we have charted operator by 
operator the relative weights of the leading order computation, of the one-loop chiral 
corrections and of the NLO 0(p 4 ) corrections. 

A 2 




Ql Q2 Q3 Q4 Q5 Q6 All 



Figure 5: Anatomy of the A2 amplitude in units of 10~ 8 GeV for central values of the input 
parameters: LO calculation (black), LO with chiral loops (half-tone), complete NLO result 
(gray). In the LO case we have taken / = /„., whereas in the remaining hystograms the 
renormalized value / = 86 MeV has been used. 



A 




Ql Q2 Q3 Q4 Q5 Q6 All 



Figure 6: Anatomy of the Aq amplitude in units of 10 GeV for central values of the input 
parameters: LO calculation (black), LO with chiral loops (half-tone), complete NLO result 
(gray). In the LO case we have taken / = f n , whereas in the remaining hystograms the 
renormalized value / = 86 MeV has been used. 

As one can see from Figs. 5 and 6, the combined effect of chiral one- loop and NLO 
xQM corrections is rather large in the case of Aq where it leads to an increase of the 
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amplitude by about 50%, while for the amplitude A 2 chiral loops are negligeable and 
xQM 0(p 4 ) corrections reduce it by 25%. It is important to remark that the final 0(p 4 ) 
result goes in the direction indicated by the AI =1/2 rule, that is, of making A larger 
and A2 smaller. 

Another remark regards the gluon penguin operators. Fig. 6 makes clear that their 
contribution to the A amplitude is not as large as often claimed. The 0(p 4 ) correction 
is almost completely accounted for by the chiral loops. 

The dependence on m c is not negligible. For instance, by varying m c (m L ) from 1.4 to 
1.3 GeV, the value of | (qq) j 1 / 3 required in order to fit the rule increases by about 10%. 

A convenient way of analizing the size of hadronic matrix elements in different the- 
oretical approaches is via the Bi factors which quantify the deviation of the hadronic 
matrix elements in a particular computation from those obtained in the VSA. According 
to the discussion below eq. (5.1), we define the Bi as 



B 



(0,2) 



Re 



(Q 



model 
«/0,2 



(Qi)o$ A 



(5.9) 



where the VSA matrix elements are by construction real. 

In Table 1 we give the Bi coefficients related to the operators Q1-6 ad, fx — 0.8 GeV. 
The three columns show how the Bi vary from LO (which includes the xQM gluonic 
corrections) to the complete NLO result, via the inclusion of meson loops. Their scale 
dependence has been already discussed in ref. [12], where also the 75-scheme dependence 
was shown. Since the NLO corrections do not affect the latter we will not repeat the 
discussion here. Anticipating the results of the next section, we have restricted the values 
of the input parameters to those required by the fit of the AI = 1/2 rule. 





LO 


+ x-l 00 P s 


+ NLO 


Bf 


4.2 


7.8 


9.5 


B?> 


1.3 


2.4 


2.9 


B? = Bf 


0.60 


0.58 


0.41 


B 3 


-0.62 


-1.9 


-2.3 


B, 


1.0 


1.6 


1.9 


B 5 ~B 6 


1.2 + 0.72 


2.0 ^ 1.2 


1.9 -=- 1.2 



Table 1: The Bi factors in the xQM including meson- loops and xQM NLO corrections. We 
have taken the gluon condensate at the central value of eq. (5.6), while the ranges given for 
-B5-6 correspond to varying the quark condensate according to eq. (5.7). The results shown are 
given in the HV scheme for M = 200 MeV and / = 86 MeV, except for the first column where 
f = f-K has been used. 

A few comments are in order. The large values for b[°1 and the small ones for Bfl 
reflect the AI = 1/2 selection rule, that is well reproduced in our approach. The decrease 
in the parameter for the operator as we increase the value of quark condensate is the 
consequence of the linear dependence on the quark condensate in the xQM with respect 
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to the quadratic one in the VSA. Finally, we confirm the rather large and negative value 
for B 3 ; the effect of such an unexpected result is discussed in ref. [12]. 

We conclude by updating in Fig. 7 the "road" to the AI = 1/2 rule already presented 
in ref. [12] were we plot the values of the amplitudes Aq and A 2 as a functions of the 
various components for the central values of the input parameters, namely M = 200 
MeV, (qq) = (-240 MeV) 3 , (a s GG/n) = (334 MeV) 4 and m s (0.8 GeV)= 220 MeV. This 
plot includes the NLO 0(p 4 ) corrections (next-to-last step) that we have here computed. 
Fig. 7 depicts in a suggestive way the decomposition of the fit into its model-dependent 
components. 



< 1.5 



(1) 

' * .(3,4) 
• • 

(5) (6) 
' • (8) 










0.5 1 1.5 2 2.5 3 3.5 

A(I=0) x 10 A 7 (GeV) 



Figure 7: The road to the AI = 1/2 rule. The numbered points are discussed in the text. 

Point (1) represents the result of free quarks (no QCD, only the operator Q 2 is present) 
combined with the VSA for the matrix element of Q 2 . Step (2) includes the effects of 
perturbative QCD renormalization on the operators Qi,2- Steps (3,4) show the effect of 
including gluon and electroweak penguin operators. The latter are responsible for the 
small shift in A 2 . Therefore, perturbative QCD brings us from (1) to (4). A crucial 
contribution is given by non-factorizable gluon condensate effects which bring us from (4) 
to (5) — still remaining at the leading 0(p 2 ) in the chiral expansion. Moving the analysis 
at the NLO, chiral loops computed on the LO chiral lagrangian lead us from (5) to (6). 
The NLO 0(p A ) corrections carachteristic of the xQM calculated in this paper yield the 
point (7). Finally, step (8) represents the inclusion of ~K-r\-v( isospin breaking effects which 
increase A 2 by about 25%. 

Let us remark that the exact match with the experimental values should not be 
taken as a theoretical prediction but rather as the proof of the reproducibility of the 
experimental rule within our approach. It should be noted once more that the bulk of 
the effect — up to point (6) — was already included in our previous analysis. 
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6 The K°-K° Mixing Parameter B K 



The scale-independent parameter Bk is defined as the product of the scale dependent 
B K ([i) parameter in 

(K \QM\K°) = p K m K B K (pi) : (6.1) 

and the function b(fi) introduced in eq. (1.7). Its determination is of crucial relevance in 
the physics of kaons. We discuss two independent ways by which it can be determined 
within the xQM approach. 

6.1 Bk from a Direct Computation 

A "model independent" estimate of Bx can be obtained by considering the relationship 
between the AS = 2 matrix element and that of the AS = 1 and AI = 3/2 amplitude 
A{K + — > 7r + 7r°) on the basis of the chiral symmetry [29] 

4 ,2 2 B ^ A m K %) A/K r+ + ox f ao\ 

- s f K <n K B K = - — - j — - A(K * 7T 7r ) . (6.2) 

Having a model that reproduces the experimental 1 = 2 amplitude, we must subtract 
from the expression of A(K + — > 7r + 7r°) all the chiral symmetry breaking components due 
to chiral loops, NLO corrections, electroweak penguins and n — r] mixing. In this way 
one obtains in the xQM approach, on the basis of chiral symmetry alone, the following 
prediction [12]: 



' : 1 ' / ! r ' (6.3) 



Bk = jb (/i) 

4 JK 



which includes the non-factorizable gluonic corrections. 

If we choose for the gluon condensate the value {a s GG/7f) = (334 MeV) 4 , which gives 
the best fit of A(K + — > 7r + 7r°) at OQo 4 ), we obtain at /i — 0.8 GeV and in the HV scheme 
the prediction 

B K ~ 0.47 , (6.4) 

to which all the specific chiral symmetry breaking corrections must be added. Notice 
that this value is about two times smaller than the LO result depicted in Fig. 8. This is 
due to the presence of f n in eq. (6.3) as derived from eq. (6.2), both in the factor f%/ f K 
and in 5{gg)- m our leading order estimate we have instead taken / = f K (f%/ fx ~ ¥ 1)> 
which we believe is the consistent choice for the tree level evaluation of the K°-K° matrix 
element (as a matter of fact it minimizes the NLO corrections). Thus the leading order 
XQM component of the 0(p 4 ) fit amounts to 

B K ~ 0.9 . (6.5) 

Adding to eq. (6.3) the meson-loop corrections (f w — > /) and the NLO corrections of 
eq. (4.19), we find the following formula 



B K = jb (ji) 

4 JK 
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B K 




LO LO+loops NLO 



Figure 8: Anatomy of Br for central values of the input parameters in the HV scheme: leading 
order result (LO), with chiral loops (LO + loops), and complete NLO result (NLO). In the LO 
case we have taken f = fx, whereas in the remaining hystograms the renormalized value / = 86 
MeV has been used. 



5(GG) 



1 + A - P (jjl) + ip r 2S {p) 



(6.6) 



where the chiral loop effects are included through the function T2s{p) of eq. (4.21). The 
functions p and A are given by 



A 



m s (p) +m(p) 



M 



6 



.M2 



mi / M z 
3M 2 ^ A 2 



(6.7) 
(6.8) 



A simpler expression of eq. (6.6) is obtained by replacing / by via the proper xQM 
and meson loop renormalizations discussed in subsections 3.2 and 3.3. We then obtain 



B K = -b(p) 



1 



1 

N~ c 



Hgg) 

N r . 



(l+\-2p(p)+r 2S (p)) 



(6.9) 



where r2s(p) is given in eq. (4.22). This corresponds to eq. (3.17) in ref. [17] with the 
addition of the 0(p A ) xQM corrections here computed. 

By varying (a s GG/Tt) in the range of eq. (5.6), M in the range of eq. (5.5) and Aq^ d 
in that given in eq. (1.4), we find 



B K = 1.1 ±0.2, 



(6.10) 



where the error includes a 8% uncertainty due to the 75-scheme dependence of b(p) at 
p = 0.8 GeV. 

In Fig. 8 we have charted the relative weight of the LO, chiral loop corrections and 
NLO contributions for central values of the input parameters.. The renormalization of Br 
from 0.9 to 1.1 amounts to about 20% and it is due to the effect of the NLO corrections 
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336 <GG> 

0.338 



Figure 9: The scale independent parameter Bk as a function of M and (GG) =(a s GG /V) 1 / 4 
(in GeV) for m<j(0.8 GeV) = 200 (lower surface) and 240 MeV (upper surface). 

that increase the value of the parameter (having chosen / = fx in the LO estimate 
minimizes the 0(p A ) corrections). This pattern differs from that found in ref. [28] where 
the NLO corrections act in the opposite direction, reducing the final value. However 
one has to keep in mind that the anatomy of the various contributions depends on the 
subtraction prescription used in the renormalization of the meson loops. 

To give a direct impression of the variation of Bk as we vary the input values, we 
show in Fig. 9 the overall dependence on (a s GG/7r) , M and m s . We recall that the NLO 
result for Bk is very little dependent on the value of / used, in the range of eq. (3.12), 
provided a corresponding fit of the A J =1/2 rule is performed. 

This computation has a systematic scale uncertainty — already present at the leading 
order, and discussed in ref. [17] — due to the poor matching between the long- and the 
short-distance computations, at variance with the cases of the amplitude A and e'/e [13], 
where the scale dependence of the Wilson coefficients is drastically reduced by the meson 
loops. This is a consequence of the chiral-loop scale dependence of as2 (/i) that adds 
coherently to the running of &(//). In this respect higher order corrections do help, since 
the running of m s improves the scale stability. Changing the matching scale from 0.8 to 
1 GeV has the effect of increasing the value of B K by about 15%. 

As a final remark, in the chiral limit where m q = (but vtik ^ 0), the NLO corrections 
change sign and actually decrease the final value of the parameter; in this case, we find 

B K (m q = 0) = 0.47 ±0.11. (6.11) 

This dramatic change qualitatively agrees with what found in ref. [30]. 

6.2 B K from the K L -K S Mass Difference 

An alternative and independent way of deriving Bk is by computing the K^-Ks mass 
difference 

AM LS = m L ~m S} (6.12) 
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and then comparing with the experimental value 

AMf s p = (3.510 ±0.018) x 10~ 15 GeV. (6.13) 

This computation is made of two parts: 

AM LS = AM SD + AM LD , (6. 14) 

where AM S d accounts for the short- distance contribution due to the AS = 2 effective 
lagrangian in eq. (1.5) and AM LD for the long-distance effects generated by the double 
insertion of the AS = 1 chiral lagrangian [17]. 

By combining eq. (6.13) and eq. (6.14) we can extract the value of B k as 

x AMff — AM LD 

8 "' M w< (6 - 15 » 

where Cs 2 (/-0 is defined in eq. (1.5). 

In a previous paper [17] we have calculated within the xQM the long-distance con- 
tributions to AMls and shown that the net result has the sign opposite to that of the 
dominant short-distance component, so that the long-distance corrections to the mass 
difference actually make the prediction smaller 1 . 

According to eq. (6.15), the value of B K that fits the experimental value of AM LS is 
found to be 

B K = 1.2 ±0.1, (6.16) 

where the error is obtained by a flat span of the whole range of the input values (Aq^d 
included). The effect of changing the matching scale from 0.8 GeV to 1 GeV reduces the 
absolute value of AM LD) corresponding to a decrease of B K by less than 10%, which is 
within the error given. 

The overlap of the two results (6.10) and (6.16) is rather encouraging and should 
dispell the concern about the rather large value of B K we find. 
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A 0{m 2 q ) Contributions to (Qq) 

As pointed out in ref. [14] the matrix element of Qq can be written as 

(7r + 7r-|Q 6 |ir°> = 2(7r-\u l5 d\0}(7v + \su\K } -2 (n + n-\dd\0}(0\s-f 5 d\K°) 

±2 [(0|ss|0) - (0|dd|0)l {7i + 7i-\s^d\K } . (A.l) 



1 In the present estimate of AMld we have not used the Gell-Mann-Okubo (GMO) relation (as we 
did in ref. [17]). As a consequence AM^jj includes also the single-pole diagrams that vanish in the limit 
in which the GMO relation holds. 
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The last line in eq. (A.l) was previously neglected because (-k + tt \s^ 5 d\K°) is zero to 
second order in momentum. However, there are 0{m 2 q ) terms, obtained when [(0|ss|0) — 



(0|dd|0) 



(m s — ma) is combined with 0{m q ) contributions from (tt + tt \s , ~f 5 d\K ), 



which a priori have to be included in a NLO analysis. A direct calculation gives 



{ ^-\s l5 d\K°) = -iV2^ + f -± 



(m s + 7n d ) + (m d + 7n u ) 



(A.2) 



and we must worry about terms proportional to m s that can be numerically relevant. 
However, some parts of this contribution can be rotated away, as we shall explain below. 
In fact, only the last term proportional to m d + m u corresponds to a physical contribution. 

The quark condensate term in eq. (A.2) combined with the m s — m d component of 
(0| ss 1 0) — (0\dd\0) correspond to a chiral lagrangian term proportional to 



Tr 



YtM a \- + \-MtV 



(A.3) 



where A± = (A6 ± M)/2 are the combination of Gell-Mann matrices projecting out AS = 
±1 transitions. This term contains K° to vacuum transitions and can therefore be rotated 
away [14] in agreement with the FKW theorem [31]. In general, 



(0\Q 6 \K° 



= 2 



(0|ss|0) - (0|dd|0)l (0\s-f 5 d\K° 



(A.4) 



and we have to find out whether there are also terms 0(m q ) which can be rotated away. 
To investigate this we do a calculation in the rotated picture of the xQM. Within this 
picture [14], 



Q 6 = -&{F { _ ) ) af} {Q L ) a {Q R ) s {Q R ) s {Q L ) fi , 



(A.5) 



where F(_) = £A_£ T , and the Greek letters are flavor indices for the constituent quark 
fields Ql = iq L and Q R = gq R . Here £f = E, L = (1 - 75 )/2, and R = (1 + 75 )/2. To 
obtain the contribution O(m^) from Q 6 to the chiral lagrangian we have to contract the 
quark fields in Q 6 to the product of two quark loops with two mass insertions. Within 
the rotated picture, the mass term reads C mass = — QM q Q, where 



M q = M q £} L + ZMttR 



Mg + 75 



(A.6) 



where and are defined in an obvious manner to be independent of 75, and A4 q 
is the current mass matrix. There are three diagrams. First, there is a diagram (A) with 
two mass insertions in the loop involving the right-handed vertex, and no mass insertion 
in the second quark loop (which corresponds to the quark condensate divided by two), 
giving the contribution : 



(qq) (-iN c ) Tr F { ^RS(p) M q S(p) M q S(p) 



(A.7) 



where S(p) = (7 • p — M) 1 , and the trace is both in flavor and Dirac spaces. Further 
there is a diagram (B) with both mass insertions in the loop with the left-handed vertex, 
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and no mass insertion at the right-handed vertex, and finally a diagram (C) with one 
mass insertion in both loops. The sum of the chiral lagrangians for the diagrams A,B,C 
are a linear combination of the terms 



Cxy = 4 Tr 



F{-\ M? M Y 



with (XY) = (SS), (PP), (SP), (PS). We have found that the term 

Cpp = Tr [A_7WjSA^;s + \-YtM q T$M q - \-YtM q M\Y. 



(A.8) 



(A.9) 



contains no K to vacuum transitions and will correspond to physical effects. In the total 
lagrangian 0(m 2 ) we also have terms involving the combinations 



C 



ss 



Cpp -\- ISC , £*sp — C 



PS 



2SC , 



where 



SC = Tr 



A_ M a Ml li 



(A.10) 



(A.11) 



This term has a K° to vacuum transition, and can be rotated away. Therefore, we obtain 
the total contribution 0(m 2 ) from Q Q as 



C-a + Cp + C-c — 



(gg) 
M 



+ 



M 2 



f 2 C PP + G S 5C , 



(A.12) 



where the term proportional to 5C has no physical consequences. 

What we found is that the terms in eq. (A. 2) containing the numerical factor 1/3 
correspond to SC. Thus, the numerically most important terms ~ m 2 are rotated away. 
Adding the last term proportional to m u + in (A. 2) to what we obtain from the 
building blocks in section 2, we obtain the result contained in the quantity 7 in (4.15), 
which agrees with the result in eq. (A.12). 



B Computing L5 and Lg in the Strong Sector 

The higher order (renormalized) counterterms in the strong chiral lagrangian can be 
determined by using the xQM. This has been done in [10] using the technique of path 
integral; here we re-derive some of them in a different and simpler way. 

The method consists in computing the Green functions with mesons as external states 
by using the xQM lagrangian given in eq. (1.9) and then to compare it with the ones 
computed using the strong chiral lagrangian 2 . From this comparison the higher order 
counter terms appearing in the strong chiral lagrangian can be determined. This is also 
the way we have followed in our previous work [11] to determine the AS = 1 weak chiral 
lagrangian up to 0(p 2 ). 

2 It should be understood that the computation includes only tree diagrams induced by 0(p 2 ) and 
0(p 4 ) terms in the strong chiral lagrangian, since the chiral loop contributions are common for the two 
approaches. 
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As an example, we present the determination of L 5 and L 8 . The same method can be 
extended to determine the other higher counter terms in the strong chiral lagrangian. 
The propagator of the meson field can be represented in momentum space as, 

(B.l) 



p 2 — Wq — E (p 2 ) ' 



where E (p 2 ) is the self-energy of the meson and m is its tree- level (bare) mass. The 
self energy E (p 2 ) may be organized as an expansion around an arbitrary point in the 
momentum space //, 



E (p 2 ) = E (/?) + (p 2 - tf) E' (/i 2 ) + E (p 2 ) , 



(B.2) 



where the prime indicates the derivative with respect to p 2 and E (p 2 ) denotes the other 
terms in the expansion. Fixing /i to be the on-shell momentum, and defining the physical 
mass as the pole of the propagator, eq. (B.l) and eq. (B.2) lead to 



ml + E (m 2 ) = m 2 , 



(B.3) 



which relates the physical mass and the bare mass. The wave-function renormalization 
is then given by 



Z = 



1 - E' ( m 



E' [m 



+ 



= 1 + E' (m 2 

Computing the two-point function by means of the xQM, we get 



(B.4) 



M(m u + m d ) 2 2 
-p - (m u + m d ) 



P 

A 2 

x 



A 2 

x 



(m u + m d ) 



m u m d 2 



10M 2 A 2 



: r "> ! - 6^=p^ + O (mj) 



MA 2 



A 2 

x 



(B.5) 



From eq. (B.3) and eq. (B.4), we obtain 



^ = 1 "2^ + 6 



M (m u + m d ) m u m d (m u + m d ) m 2 



and 



A x 



+ 2- 



3mf 



MA 2 



10M 2 A 2 



+ 



A 2 
. x 



M (m u + m d ) 



A 2 

x 



- 2 



(B.6) 



with 



0\2 



< = K) 



(m u + m d y 

A 2 

x 



5 2 + 6 



M (m M + m d ) z 

A 2 

x 



5 - (m u + m d y 



(m u + m d 



MA 2 



5o 



fi 2 - 6 ^(^ + ^) Bo + o « 



rn, 



p 



mz 



(m u + m d ) 



(B.7) 



(B.8) 
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The decay constants f n and fx are going to be necessary in determining L 5 . The coupling 
f n computed in the chiral quark model is given at 0(p 2 ,m q ) by 



U = f 



(qq) 2M ' 



^x 



(B.9) 



The results for f n and m\ computed from the strong chiral lagrangian by including the 
relevant 0{p A ) counterterms (but no chiral loops) are, 



U = f 



A mi. T 8m 

~F 5 + — P 



— Li A 



and 



mi = {mlY + -j^Bl (m u + m d ) z [2L 8 - L 5 ] 



Comparing eq. (B.9) and eq. (B.10) implies 

5 8(qq) M \ 



1-6 



(B.10) 



(B.H) 



(B.12) 



which is also obtained by a direct calculation within the xQM (in the rotated picture of 
[14]). Accordingly, at the same level of approximation, we find 



L 4 = 0. 



From eq. (B.7) and eq. (B.ll) we get at 0(m 

P 



2La — Lr\ 



8BI 



B 2 M 

X X 



(B.13) 



(B-14) 



Eq. (B.14) together with eq. (B.12), yields 

N c 1 P 



U = 



16tt 2 24 lQ(qq)M 



MP 

(qq) 



(B.15) 



The coupling f K and m K could be equally used in determining L 4 , L 5 and L 8 leading 
to the same result. 
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C Input Parameters 



ngro m or pr 


~\ 1 1 1 

VcXl LLC 




n 93 


TH, 7 


91 187 GeV 


Tfl\Y 


80 33 GeV 


pole 


175 ± 6 GeV 


m i ( hi i i 
ni b\ Ui b) 


4 4 GeV 




1.4 GeV 


m s (1 GeV) 


178 ± 18 MeV 


m„ 4- t(Ia (1 GeV) 


12 ± 2.5 MeV 


A (4) 


S40 + 40 MeV 






T/ 




M 


9nn + 5 MoV 


laa\ 


— ^Z41U ivie V J 


{a s LrLl/TT) 


^Oo4 ± 4 ivie V J 


J 




fir = /tt+ 


92.4 MeV 


Ik = fx+ 


113 MeV 


m n = {m n + + m 7r o)/2 


138 MeV 


niK = m K o 


498 MeV 


m ri 


548 MeV 


cos 5q 


0.8 


COS ^2 


1.0 



Table 2: Table of the numerical values of the input parameters. 
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